A0-A178  885  MODELING  XV-15  TILT-ROTOR  AIRCRAFT  DVMAKICS  8V 

FREOUENCV  AND  TIHE-DONAIN  <U)  ARHV  AVIATION  RESEARCH 
AND  TECHNOLOGV  ACTIVITV  HOFFETT  FIELD  C 
UNCLASSIFIED  H  8  TISCHLER  ET  AL  16  OCT  86  F/G  1/2 


3 1. 


uTJC  FILE  copy 


ADVISORY  GROUP  FOR  AEROSPACE  RESEARCH  &  DEVELOPMENT 


7  RUE  AMCELLE  92200  NEUILLY  SUR  SEIME  FRANCE 


MODELING  XV- 15  TILT-ROTOR  AIRCRAFT  DYNAMICS  BY  FREQUENCY 
AND  TIME -DOMAIN  IDENTIFICATION  TECHNIQUES 


Mark  B.  Tischler 
JUrgen  Kaletka 


s 


One 

ELECTE 
APR  0  7  1987 


I 


Paper  No.  9 

AGARD  Flight  Mechanics  Panel  Symposium 

on 

ROTORCRAFT  DESIGN  FOR  OPERATIONS 
Amsterdam,  The  Netherlands 
13-16  October,  1986 


NORTH  ATLANTIC  TREATY  ORGANIZATION 


MODELING  XV- 15  TILT-ROTOR  AIRCRAFT  DYNAMICS  BY  FREQUENCY  AND 
TIME-DOMAIN  IDENTIFICATION  TECHNIQUES 

by 

Mark  B.  Tlschler 
Aerofl ight dynamics  Directorate 
US  Army  Aviation  Research  and  Technology  Activity 
Ames  Research  Center 
Moffett  Field,  California  94035-1099 

Juergen  Kaletka 
Research  Scientist 

Deutsche  Forschungs-  und  Versuchsanstalt 
fuer  Luft-  und  Raumfahrt  e.V.  (QFVLR) 

Institut  fuer  Flugmechanik 
D-3300  Braunschweig-Flughafen,  Germany 


SUMMARY 

Models  of  the  open-loop  hover  dynamics  of  the  XV-15  Tilt-Rotor  Aircraft  are  extracted  from  flight 
data  using  two  approaches:  frequency-domain  and  time-domain  identification.  Both  approaches  are  reviewed 
and  the  identification  results  are  presented  and  compared  in  detail.  The  extracted  models  compare  favor¬ 
ably,  with  the  differences  associated  mostly  with  the  Inherent  weighting  of  each  technique.  Step 
responses  are  used  to  show  that  the  predictive  capability  of  the  models  from  both  techniques  is  excel¬ 
lent.  Based  on  the  results  of  this  study,  the  relative  strengths  and  weaknesses  of  the  frequency  and 
time-domain  techniques  are  summarized,  and  a  proposal  for  a  coordinated  parameter  identification  approach 
Is  presented. 
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roll  rate,  pitch  rate,  and  yaw  rate,  respectively,  deg/sec 
longitudinal,  lateral,  and  vertical  velocities,  respectively,  m/sec 
coherence  function  between  variable  x  and  y 

aileron  surface  deflection  (deg),  elevator  surface  deflection  (deg),  and 
deflection  (deg),  respectively 

power  lever  deflection,  X 

damping  ratio 

time  delay 

roll,  pitch,  and  yaw  angles,  respectively,  rad  (deg) 
undamped  natural  frequency,  rad/sec 
Inverse  time  constant,  rad/sec 
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1.  INTRODUCTION 

Dynamics  Identification  methodologies  generally  fall  into  two  categories:  frequency-domain  and  time- 
domain.  The  choice  of  techniques  to  be  used  Is  usually  based  on  the  analyst's  personal  familiarity  with 
the  methods  and  on  the  specific  application.  Each  approach  has  Inherent  strengths  and  weaknesses. 
Frequency-domain  identification  uses  spectral  methods  to  determine  frequency  responses  between  selected 
Input  and  output  pairs.  Then,  least-squares  fitting  techniques  are  used  In  the  frequency-domain  to  obtain 
closed-form  analytical  transfer-function  models  of  linear  input-to-output  processes.  Time-domain  identi¬ 
fication  first  requires  the  selection  of  a  state-space  model  structure,  which  may  be  linear  or  nonlinear. 
Model  parameters  are  Identified  by  least-squares  fitting  of  the  response  time-histories  or  by  maximum 
likelihood  methods.  Transfer  functions  for  linear  models  and  frequency  responses  can  then  be  obtained 
from  the  identified  state-space  formulation. 


The  US  Army  has  been  developing  frequency-domain  identification  techniques  In  support  of  handling 
qualities,  flight,  and  simulation  experiments.  Extensive  flight  experiments  have  been  conducted  on  the 
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XV-15  Tilt-Rotor  Aircraft  (Fig.  1).  References  1  and  2  present  the  identified  open-loop  frequency 
responses,  transfer  functions,  and  model  verification  results.  Frequency-domain  Identification  tests  have 
also  been  recently  conducted  on  the  Bell-214-ST  single  (teetering)  rotor,  and  CH-47  tandem-rotor  aircraft 
(Refs.  3  and  4).  In  the  federal  Republic  of  Germany,  the  DFVIR  has  had  extensive  experience  with  maximum- 
likelihood,  time-domain  Identification  techniques.  Linear  and  nonlinear  model-identification  methods  have 
been  developed.  Much  of  the  DFVLR  experience  with  helicopter  identification  has  been  associated  with  the 
highly  coupled  60-105  hlngeless-rotor  helicopter  (Refs.  5  and  6). 

As  part  of  an  ongoing  US/FRG  Memorandum  of  Understanding  (MOU)  on  helicopter  flight  control,  an 
extensive  joint  study  Is  being  conducted  to  analyze  the  XV-15  data-base  for  the  (open-loop)  hover  flight 
condition  using  both  time-  and  frequency-domain  techniques.  The  primary  objectives  of  this  study  are 
to:  (1)  gain  a  better  appreciation  for  the  relative  strengths  and  weaknesses  of  each  technique;  and 
(2)  develop  Improved  methods  of  identification  for  rotorcraft. 

This  paper  reviews  the  dynamics  identification  techniques  which  have  been  developed  In  the  US  and  the 
Federal  Republic  of  Germany.  The  results  of  applying  these  techniques  to  the  XV-15  data  base  are  pre¬ 
sented  and  compared,  and  sources  of  differences  In  the  extracted  models  are  discussed.  Finally,  conclu¬ 
sions  concerning  the  appropriate  applications  for  each  technique  and  proposals  for  unified  Identification 
methods  using  both  approaches  are  presented. 


Fig.  1.  The  XV-15  Tilt-Rotor  Aircraft,  a)  Hover  configuration;  b)  cruise  configuration. 


2.  OVERVIEW  OF  FREQUENCY-DOMAIN  AND  TIME-DOMAIN  IDENTIFICATION  TECHNIQUES 

A.  Frequency-Domain  Identification  Method 

The  frequency-domain  Identification  approach  developed  by  the  US  Army  is  depicted  in  Fig.  2.  Spec¬ 
tral  methods  based  on  the  Chirp  z-transform  are  used  to  extract  high-resolution  frequency  responses 
between  selected  Input  and  output  pairs.  The  identification  results  are  presented  in  Bode-plot  format: 
magnitude  and  phase  of  the  output  to  the  input  versus  frequency.  These  identification  results  are  non- 
parametric  because  no  model  structure  has  been  assumed.  As  such,  they  can  be  very  useful  for  flight- 
control  system  design  and  handling-qualities  compliance  testing;  for  example,  currently  proposed  handling- 
qualities  criteria  for  the  LHX  (Ref.  7)  are  based  on  frequency-domain  parameters  which  can  be  read 
directly  from  these  graphical  results.  Frequency  responses  obtained  from  real-time  and  nonreal-time  simu¬ 
lations  can  be  compared  directly  with  the  flight  data  to  expose  limitations  and  discrepancies  In  the  simu¬ 
lator  models  (Ref.  1).  The  fact  that  this  comparison  can  be  made  Initially  without  an  a  priori  assumption 
of  model  structure  or  order  is  especially  Important  for  verifying  mathematical  models  of  new  aircraft 
configurations.  When  the  model  structure  and  parameter  values  are  required,  they  may  be  obtained  by  fit¬ 
ting  the  tabulated  frequency-responses  with  analytical  transfer-function  models  to  extract  modal  charac¬ 
teristics.  Examples  of  this  application  are  the  testing  of  handling-quality  specifications  gfven  In 
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Fig.  2.  Frequency-domain  Identification  method. 


lower-order  equivalent  system  terms,  and  the  examination  of  transfer  function-based  control  system 
designs.  Since  this  fitting  procedure  Is  completed  after  the  frequency  response  Is  extracted,  the  order 
of  the  transfer  function  can  be  carefully  selected  to  avoid  an  overparameterized  model.  Multi-Input/ 
multi-output  frequency-response  methods  are  suitable  for  extracting  a  transfer  matrix  which  Includes  the 
Important  coupling  effects.  Finally,  the  extracted  models  are  driven  with  the  flight-test  control  Inputs 
to  verify  the  time-domain  response  characteristics. 

The  semllog  frequency  format  of  the  Bode-plot  presentation  and  subsequent  transfer-function  fit  makes 
the  Identified  transfer-function  and  state-space  models  most  accurate  at  mid  and  high  frequency  (Initial 
time  history  transients).  The  low-frequency  and  steady-state  response  prediction  of  the  extracted  models 
Is  generally  not  as  good  as  In  time-domain  Identification  approaches. 

B.  Time-Domain  Identification  Method 

The  general  approach  used  In  time-domain  Identification  Is  shown  In  Fig.  3.  Time-based  Identifica¬ 
tion  techniques  are  Initially  applied  to  the  data  to  check  their  Internal  compatibility.  Data  inconsis¬ 
tencies  resulting  from  calibration  errors,  drifts,  or  Instrumentation  failures  are  detected  by  comparing 
redundant  measurements  from  Independent  sensors,  such  as  rate  and  attitude  gyros,  or  altitude  change  and 
vertical  acceleration  (Ref.  6).  This  approach,  which  can  be  used  on-line,  helps  to  ensure  that  only  con¬ 
sistent  data  are  generated  for  the  further  evaluation  and  system  Identlf Icatlon. 

For  this  next  step,  the  aircraft  dynamics  are  modeled  by  a  set  of  differential  equations  describing 
the  external  forces  and  moments  In  terms  of  accelerations  and  state  and  control  variables.  The  coeffi¬ 
cients  In  these  equations  are  the  stability-derivatives.  In  some  cases,  a  priori  values  for  these  deriva¬ 
tives  can  be  obtained  from  analytical  calculations,  wind-tunnel  data,  or  from  start-up  Identification 
techniques  such  as  a  least-squares  method.  The  responses  of  the  model  and  aircraft  resulting  from  the 
flight-test  control  inputs  are  then  compared.  The  response  differences  are  minimized  by  the  identifica¬ 
tion  algorithm  which  iteratively  adjusts  the  model  parameters.  In  this  sense,  aircraft  system  identifica¬ 
tion  Implies  the  extraction  of  physically  defined  aerodynamic  and  flight  mechanics  parameters  from  flight- 
test  data.  Usually,  it  Is  an  off-line  procedure  since  some  skill  and  Iteration  are  needed  to  develop  an 
appropriate  model  formulation.  Model  formulation  Involves  consideration  of  model  structure,  selection  of 
significant  parameters,  and  Inclusion  of  Important  nonllnearltles.  Tfme-domaln  techniques  yield  a  multl- 
Input/mu It  1 -output  model  that  Is  appropriate  for  application  In  stability  and  control  analysis,  simula¬ 
tion,  and  control  system  design.  The  Identified  parameters  are  also  useful  for  comparison  and  correction 
of  analytically  or  wind-tunnel  derived  stability-derivatives  (Ref.  8). 
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Fig.  3.  Time-domain  identification  method. 


A  key  feature  of  the  time-domain  identification  technique  Is  that  the  extracted  models  are  based  on 
the  curve  fitting  of  the  original  measured  (time-domain)  flight-test  data.  Errors  which  may  occur  In  the 
transformation  of  the  data  from  the  time-  to  the  frequency-domain  are  thus  avoided.  The  Identified  models 
can  then  be  easily  presented  In  the  frequency-domain  as  Bode  plots  or,  if  the  Identified  model  Is  linear, 
also  as  parametric  transfer-functions. 
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3.  IDENTIFICATION  OF  XV-15  OPEN-LOOP  DYNAMICS  IN  HOVERING  FLIGHT 

This  section  reviews  the  XV-15  flight-test  data  base  for  parameter  Identification  and  presents  and 
compares  the  results  of  (linear)  frequency  and  time-domain  identification  methods  for  the  open-loop  hover 
flight  condition.  For  Illustrative  purposes,  the  roll  response  Identification  Is  discussed  In  detail. 


A.  Flight-Test  Data  Base 

The  complete  data  base  for  dynamics  Identification  Includes  four  flight  conditions  from  hover  to 
cruise.  The  present  study  concentrated  exclusively  on  the  identification  of  the  open-loop  hover  dynamics 
because: 

The  dynamics  for  this  flight  condition  are  coupled  and  very  unstable  which  makes  this  case  the  most 
difficult  to  analyze. 

Nonlinear  effects  are  the  most  significant  in  the  hover  flight  condition,  which  allows  a  good  demon¬ 
stration  of  the  nonlinear  Identification  techniques  developed  by  the  QFVLR. 

Focusing  on  the  rotor-borne  flight  condition  maximizes  the  carry-over  of  the  present  experience  to 
future  rotorcraft  identification  studies  to  be  carried  out  under  the  MOD. 

The  pitch  and  roll  axis  Instabilities  for  the  hover  flight  condition  are  characterized  by  a  time-to- 
double  amplitude  of  about  3  sec.  Therefore,  long-period  inputs  needed  to  identify  the  low-frequency 
vehicle  dynamics  are  not  practical  for  the  open- loop  hovering  vehicle.  Extraction  of  the  open-loop 
vehicle  dynamics  from  closed-loop  testing  is  possible  subject  to  an  Important  condition:  the  total 
surface  deflection,  which  Is  comprised  of  Inputs  from  the  pilot  and  the  stability  and  control  augmentation 
system  (SCAS),  must  contain  a  significant  component  which  Is  uncorrelated  with  the  response  of  the  vehicle 
(Ref.  9).  Then,  the  required  low-frequency  Inputs  can  be  conducted  on  the  closed-loop  (stable)  vehicle. 

Flight-Test  Inputs.  Two  types  of  Inputs  were  executed  in  the  identification  flight  tests. 
'Frequency-sweep"  Inputs  were  used  for  model  extraction,  and  step  Inputs  were  used  for  model  verification. 

Two  typical  concatenated  lateral  stick  frequency-sweeps  completed  during  the  hover  flight  tests  of 
the  XV  15  are  shown  in  Fig.  4a.  These  tests  used  pilot-generated  rather  than  computer-generated  Inputs. 
The  sweep  Is  Initiated  with  two  low-frequency  Input  cycles  corresponding  to  the  lower  bound  of  the  fre¬ 
quency  range  of  primary  Interest  (0. 2-6.0  rad/sec).  These  cycles  ensure  good  excitation  of  the  low- 
frequency  vehicle  dynamics.  After  the  Initial  two  low-frequency  cycles,  the  lateral  stick  Is  oscillated 
at  progressively  higher  frequencies  for  an  additional  50  sec.  By  the  end  of  the  90  sec  duration  test,  the 
stick  is  being  driven  at  fairly  high  frequencies  (4  Hz  shown  in  Fig.  4).  The  input  amplitudes  are  fairly 
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Fig.  4.  Two  lateral  stick  frequency-sweeps  In  hover,  a)  Lateral  stick  Inputs;  b)  aileron  surface 

deflections. 
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Fig.  5.  Roll-rate  response  (p)  during  lateral  frequency-sweeps. 

small  at  low  frequency  where  vehicle  motions  are  considerable,  with  larger  Inputs  at  mid  frequency,  and 
smaller  Inputs  again  at  very  high  frequency.  The  associated  aileron  surface  deflection  (total  Input  to 
the  aircraft*)  shown  In  Fig.  4b  reflects  a  significant  component  from  the  pilot  Input  (note  that  s.  and 
are  defined  with  opposing  sign  conventions).  The  resulting  roll-rate  amplitudes  of  10-20  deg/sec  as 
shown  In  Fig.  5  are  typical  for  frequency- sweep  tests.  The  frequency-sweep  Is  especially  well  suited  for 
frequency-domain  Identification  because  It  Is  a  periodic  input  form  that  excites  the  vehicle  in  ail  of  its 
dominant  modes  of  motion  within  the  frequency  range  of  Interest.  This  Input  also  has  some  advantages  for 
time-domain  Identification.  Vehicle  excitation  Is  restricted  to  be  within  the  frequency-range  of  model 
applicability  which  Is  especially  important  for  meaningful  state-space  parameter  results  (Ref.  4).  Also, 
the  monotonic  Increase  In  frequency  allows  the  time-domain  Identification  to  be  frequency-weighted  which 
compensates  for  the  Inherent  low-frequency  weighting  of  this  method. 

Step  Inputs  are  commonly  used  In  the  flight  test  coanunlty  to  expose  dominant  vehicle  characteris¬ 
tics,  so  they  represent  a  good  test  of  the  Identified  model's  predictive  capability.  Step  Inpjts  were 
executed  In  both  the  open-  and  closed-loop  condition.  Open-loop  verification  ensures  that  the  identified 
models  reflect  the  dynamics  of  the  open-loop  vehicle  and  not  those  of  the  Inverse  feedback  element 
(Ref.  9).  Step  Inputs  with  the  fl ight -control  system  engaged  are  also  useful  since  the  steadier  Initial 
conditions  allow  fine  differences  between  the  model  and  the  flight  responses  to  be  exposed. 

B.  Frequency-Domain  Identification 

The  most  Important  step  In  the  frequency-domain  Identification  procedure  Is  the  extraction  of  accu¬ 
rate,  high-resolution  frequency  responses  between  the  various  Input  and  output  pairs.  A  key  metric  for 
assessing  the  quality  of  the  frequency-response  Identification  Is  the  coherence  function  (TJ  ).  This 
frequency-dependent  function  Indicates  that  fraction  of  the  output  response  which  Is  linearly  related  to 
the  excitation  signal.  The  random  error  associated  with  the  frequency-response  Identification  Is  depen¬ 
dent  on  the  value  of  the  coherence  function  at  each  frequency,  and  on  the  number  of  (Independent)  time 
history  segments  ("windows,"  Nd): 

^Although  the  aileron,  elevator,  and  rudder  surfaces  are  not  effective  In  hover,  they  continue  to  be 
actuated  In  addition  to  the  primary  effectors  which  are  the  rotor  collective  and  swashplate  controls.  It 
was  found  to  be  most  expedient  to  refer  all  the  transfer  functions  to  these  surface  deflections,  since 
neglecting  the  small  servo  lags,  these  are  related  to  the  sum  of  the  pilot  and  SCAS  inputs  through  a 
mixing  ratio  which  Is  constant  across  the  entire  flight  envelope. 
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The  length  of  the  window  (secs)  determines  the  amount  of  low-frequency  power  and  the  associated  low- 
frequency  coherence  which  can  be  achieved.  Low  variance  in  the  spectral  Identification  therefore  requires 
high  coherence  and  multiple  concatenated  time-history  records.  Initial  analyses  of  the  XV-15  data  base 
used  all  available  repeat  runs  (three  were  used  in  the  original  analysis  of  Ref.  1),  without  concern  of 
the  Individual  coherence  quality  of  each  run.  Subsequent  time-domain  analyses  by  the  DFVLR  and  frequency- 
domain  analyses  by  the  US  Army  indicated  that  some  of  the  frequency-sweep  runs  were  unsuitable  for  Identi¬ 
fication,  and  should  be  removed  from  the  concatenation  procedure.  In  the  case  of  the  lateral  axis  fre¬ 
quency-sweep,  one  of  the  three  runs  was  found  to  be  unsuitable  because  of  low  coherence.  The  frequency- 
response  obtained  with  the  remaining  two  (good)  runs  has  substantially  Improved  spectral  quality.  This 
frequency  response  and  the  associated  coherence  function  are  shown  In  Figs.  6  and  7.  Good  Identification 
Is  achieved  over  the  frequency  range  0. 2-9.0  rad/sec. 

The  magnitude  response  peak  is  due  to  the  dominant  roll  modes  which  are  In  the  frequency-range  of 
0.5-1  rad/sec;  the  associated  phase  rise  Indicates  that  the  modes  are  unstable.  At  the  higher  frequencies 
(1.0-10  rad/sec),  the  magnitude  and  phase  plots  follow  a  K/s  characteristic.  The  value  of  the  constant 
(K)  Is  the  roll  response  sensitivity  (L4a).  The  relatively  flat  phase  response  at  high  frequency  indi¬ 
cates  a  very  small  value  of  effective  time  delay.  Finally,  the  drop  in  coherence  function  near  the  magni¬ 
tude  peak  suggests  the  existence  of  nonl Inearltles  for  large  vehicle  motions. 

(1)  Lateral/Olrectlonal  Transfer-Function  Models 

The  selection  of  the  order  and  structure  of  the  transfer-function  models  Is  predominantly  based  on 
three  important  factors  (Ref.  4): 

(a)  The  models  must  be  appropriate  to  the  frequency-range  of  concern  (0. 2-6.0  rad/sec  In  the 
present  study) . 

(b)  The  models  must  provide  a  reasonable  fit  of  the  Input-to-output  frequency  response  within 
the  frequency  range  associated  with  good  coherence. 

(c)  The  selected  models  should  be  based  on  a  theoretical  analysis  of  the  effective  physical 
order  of  the  system.  Therefore,  the  appropriate  transfer-function  models  are  a  function  of  flight  condi¬ 
tion  and  flight-control  system  status  { 1 .e . ,  SCAS-on  or  SCAS-Off). 

For  the  open-loop  XV-15  In  the  hovering  flight  condition,  the  yaw  (and  heave)  responses  are  essen¬ 
tially  decoupled  and  first  order  In  nature.  Therefore,  an  appropriate  model  for  yaw-rate  response  to 
pedal  Inputs  Is: 

-i,s 
"a  *  * 

T  (s)  =  7TVT7TT  (?) 


The  un-axis  roll  (and  pitch)  responses  are  dominated  by  the  hovering  cubic,  and  as  seen  In  Fig.  6  have  one 
excess  pole  at  high  frequency: 


(+> 


The  dominant  source  of  coupling  In  the  open-loop  configuration  is  the  yaw  response  to  lateral  stick 
inputs.  This  coupling  arises  from  the  rotor  torque  differential  which  accompanies  the  differential 
collective  inputs  used  for  roll  control,  frequency-response  Identification  of  the  coupled  response 
(Ref.  1)  Indicates  an  appropriate  transfer-function  model  of: 

'V 

N6  (‘/VlV“*le 

ra  <s>  1  T?/Ty)(i7Tr)ip-;;vrT"  (4) 

The  denominator  parameters  of  the  lateral/directional  transfer-function  modeis  (Eqns.  2-4)  represent 
natural  dynamics  modes  of  the  vehicle.  Therefore,  the  common  modes  must  have  the  same  values  for  all 
three  responses.  Maintaining  this  relationship  Is  essential  for  achieving  unique  and  physically  meaning¬ 
ful  transfer-function  models.  While  it  Is  possible  to  fit  all  three  responses  simultaneously  to  maintain 
the  conmonallty  of  denominator  parameters,  this  approach  Is  not  the  best.  A  better  strategy  is  to  iden¬ 
tify  individual  parameters  from  the  on-axis  frequency-response  in  which  they  have  the  dominant  effect. 

♦Window  overlapping  further  reduces  the.random  errop  below  that  shown  In  Eqn.  (1)  (Ref.  9). 

‘Shorthand  notation:  |  c ,  u)  Implies  s'  *  2;wS  +  a‘,  r.  -  damping  ratio,  -  undamped  natural 


frequency  (rad/sec);  and  (1/T)  implies  s  ♦  (!/'i,  rad/sec. 


undamped  natural 
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FREQUENCY  DOMAIN  IDENTIFICATION 


FREQUENCY,  rad/sac 

Fig.  6.  Frequency-domain  identification  of  roll- 
rate  response  to  aileron  (p/-«a).  a)  Magnitude; 
b)  phase. 


FREQUENCY,  rad/sac 

Fig.  7.  Coherence  function  for  roll-rate  response 
Identification  (r^  p). 


Then  these  parameters  are  fixed  In  the  Identification  of  the  off-axis  transfer-functions.  So,  for  exam¬ 
ple,  the  following  yaw  response  transfer-function  Is  obtained  from  the  pedal  sweeps: 


7-  <S>  = 


0.619  e 


-0.0210s 


(s  +  0.102) 


This  result  shows  that  the  yaw  response  of  the  tilt-rotor  configuration  Is  very  lightly  damped,  as  com¬ 
pared  to  a  standard  helicopter  with  a  tall  rotor.  The  small  effective  time  delay  indicates  that  lags 
caused  by  unmodeled  high-frequency  dynamics  are  negligible. 

The  next  step  Is  to  Identify  the  roll-rate  transfer  function  (p/6a).  Since  the  yaw  mode  has  been 
Identified  In  Eqn.  5,  this  parameter  Is  fixed  In  the  roll-response  transfer  function  (Eqn.  3).  Then  the 
remaining  parameters  are  varied  to  obtain  the  best  least-squares  fit  of  the  roll  rate  frequency-response 
(Fig.  6): 


The  fitting  range  Is  from  0. 2-9.0  rad/sec,  In  which  the  coherence  function  (y‘  p)  Indicates  good  spectral 
accuracy.  As  expected  from  the  phase  response  characteristics  (Fig.  6b),  the  open-loop  roll-response 
dynamics  are  dominated  by  an  unstable  roll  mode  with  the  frequency  of  about  0.4  rad/sec.  The  associated 
tlme-to-double  amplitude  Is  3.5  sec.  The  pole-zero  pair  ( 1/T„ ^) /( 1/Ty)  Is  at  very  low  frequency  and 
nearly  cancels  out.  This  reveals  that  yaw  coupling  does  not  noticeably  affect  the  roll-response  charac¬ 
teristics.  Therefore,  a  lower-order  roll  response  model  which  contains  only  the  hovering  (lateral)  cubic 
roots  (1/Tr) ltr,wr)  and  entirely  Ignores  yaw  coupling  Is  an  appropriate  approximation  for  this  vehicle. 
This  assumption  Is  common  for  hovering  aircraft.  The  low-frequency  numerator  factor  (1/Ta.)  associated 
with  lateral  translation  damping,  Is  marginally  unstable  (t'me-to-double  amplitude  =  7.5  sec)  Indicating  a 
very  low  value  of  the  velocity  damping  derivative  (Yy).  Finally,  the  effective  time  delay  for  the  roll 
response  (O  Is  small,  suggesting  that,  as  In  the  yaw  response  (Eqn.  5),  the  unmodeled  high-frequency 
lags  are  not  significant. 

With  the  lateral/dlrectlonal  denominator  factors  (dominant  vehicle  modes)  Identified  using  the 
on-axis  frequency-responses,  the  numerator  factors  of  the  off-axis  response  (r/6  )  can  now  be  extracted. 
The  denominator  factors  of  Eqn.  4  are  fixed  and  the  least-squares  fit  gives: 


In  the  frequency  range  u  >  1.2  rad/sec,  the  yaw  response  to  aileron  Inputs  Is  dominated  entirely  by  the 
coupling  derivative,  N4a.  At  low  frequencies,  the  dynamics  are  affected  by  the  unstable  hovering  cubic. 

As  shown  In  Figs.  6a  and  6b,  the  transfer-function  model  of  Eqn.  6  Is  a  good  representation  of  the 
Identified  roll  response  In  the  range  of  satisfactory  coherence  (0. 2-9.0  rad/sec).  Although  the  present 
transfer-function  model  (Eqn.  6)  Is  not  significantly  different  from  that  obtained  previously  (Ref.  1) 
using  all  of  the  available  sweep  runs  (Including  the  poor  quality  runs),  the  match  between  the  model  and 
flight  data  Is  significantly  Improved. 
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Close  examination  of  fig.  6  shows  that  the  match  between  the  model  and  flight  data  is  much  better  in 
magnitude  than  in  phase.  This  Is  because  of  the  relative  weighting  selected  for  the  identification  (1  OB 
magnitude  error:  7°  phase  error)  which  Is  common  for  lower-order  equivalent  system  matching  (Ref.  10). 

On  the  basis  of  the  steep  phase  response  of  the  flight  data  at  the  dominant  mode  (u  =  0.5  rad/sec)  as 
compared  to  the  transfer-function  model,  a  lower  (less  negative)  damping  ratio  is  Indicated.  This  incon¬ 
sistency  between  the  magnitude  and  phase  responses  indicates  nonlinear  behavior  in  the  dominant  modes  of 
roll  motion.  As  mentioned  previously,  this  is  also  reflected  by  the  drop  In  coherence  In  the  same  fre¬ 
quency  range.  Significant  side-by-side  nonlinear  rotor  interactions  are  known  to  exist  for  large  lateral- 
velocity  transients,  as  were  encountered  during  the  low-frequency  Inputs. 

The  lateral/directional  transfer-function  model  results  of  this  section  are  summarized  in  Table  1. 

(2)  Longitudinal  Transfer-function  Models 

In  hovering  vehicles,  pitch  and  roll  dynamics  are  analogous.  The  pitch  response  is  dominated  by  a 
longitudinal  hovering  cubic  analogous  to  the  lateral  hovering  cubic,  and  a  first-order  heave  response  is 
analogous  to  the  first-order  yaw  response.  Power  lever  (vertical  control)  and  longitudinal  stick  Inputs 
do  not  Induce  significant  inter-axis  coupling  in  the  XV-15  configuration. 

Spectral  analysis  of  the  individual  pitch-sweeps  showed  that  only  one  of  the  three  runs  had  satisfac¬ 
tory  coherence  for  use  in  identification.  (The  original  analysis  of  Ref.  1  used  all  three  runs.) 
Transfer-function  models  are  extracted  from  the  identified  frequency  responses  using  the  same  approach 
discussed  above  for  the  lateral/directional  dynamics.  The  heave  response  is  determined  first  from  collec¬ 
tive  sweeps,  and  then  the  pitch  response  is  determined  with  the  identified  heave  mode  ( 1 /Th)  fixed.  The 
resulting  transfer  functions  for  pitch  rate  and  vertical  acceleration  responses  are  summarized  in  Table  1. 
As  in  the  roll  case,  the  pitch  response  is  dominated  by  a  hovering  (longitudinal)  cubic,  comprised  of  a 
low-frequency  unstable  oscillation,  Up,un),  and  a  stable  aperiodic  mode  (1/Tp).  Also,  the  pitch 
transfer-function  model  fits  the  Identified  frequency  response  much  better  in  magnitude  than  in  phase. 
Based  on  phase  response  considerations  alone,  the  unstable  damping  ratio  would,  as  before,  be  much  lower 
(less  negative).  The  discrepancy  between  the  magnitude  and  phase  fit  is  again  due  to  nonlinearities 
associated  with  the  large  velocity  perturbations  encountered  during  the  low-frequency  Inputs. 

C.  Time-Domain  Identification 

Maximum  likelihood  (ML)  technique  is  generally  accepted  as  one  of  the  most  suitable  time-based 
methods  for  aircraft  parameter  Identification.  The  main  advantages  of  the  ML  estimation  are: 

(I)  It  yields  asymptotically  unbiased  and  consistent  estimates  for  linear  systems. 


TABLE  1  Comparison  of  Transfer-Function  Models  for  Hover 
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(2)  It  provides  the  Cramer -Rao-Bound,  which  is  a  measure  of  the  reliability  of  each  estimate. 

(3)  It  yields  the  correlation  between  the  identified  parameters. 

8oth  the  Cramer-Rao-8ound  and  the  parameter  correlation  help  to  develop  an  appropriate  model  struc¬ 
ture  and  to  avoid  "over-parameterization."  A  nonlinear  maximum  likelihood  method  developed  by  OFVLR 
(Refs.  11,  12)  was  utilized  for  time-domain  identification  of  the  XV - 1 5-  This  technique  allows  a  general 
linear  and  nonlinear  model  formulation  of  the  state  and  measurement  equations: 


x(t)  =  f(x(t),  u(t)  -  au,  a) 
y(t)  =  g(x(t),  u(t)  -  au,  a)  +  Ay 


x(t  =  0)  =  Xn 


x  -  computed  state  vector 
y  =  measured  variables 
u  =  measured  control  vector 
xg  =  initial  conditions 
a  =  system  parameters 
au ,Ay  =  zero  shifts 

The  initial  conditions  and  zero  shift  terms  are  included  to  compensate  for  drifts  and  offsets  in  the 
measurements. 

For  the  XV- 15  data  evaluation,  the  ML  program  was  first  utilized  to  check  compatibility  of  the  mea¬ 
surements,  and  to  reconstruct  the  nonmeasured  data.  Then  the  program  was  used  for  the  parameter  identifi¬ 
cation  Itself.  In  the  following  sections,  these  steps  are  discussed  in  detail. 

a.  Data  Compatibility  and  Reconstruction 

The  XV-15  instrumentation  system  provides  attitude  rates,  attitude  angles,  and  linear  accelerations; 
speed  measurements  for  the  hover  flight  condition  are  not  available.  Therefore,  only  the  compatibility  of 
the  angular  data  could  be  evaluated.  A  satisfactory  agreement  between  calculated  and  measured  angles  was 
found  and  no  additional  corrections  were  made.  For  the  frequency  sweeps,  speed  components  were  derived  by 
integrating  the  measured  linear  accelerations.  Since,  for  these  tests,  the  aircraft  is  In  trim  at  the 
beginning  and  the  end  of  each  sweep,  speed  equation  biases  can  be  estimated  to  meet  the  boundary  condi¬ 
tions:  u(0)  »  v(0)  =  w(0)  =  u(tp)  =  v(tp)  =  w(tp)  =  0.  For  the  system  identification,  the  calculated 
velocity  variables  are  Included  In  the  measurement  vector  together  with  the  linear  accelerations. 

Strictly  speaking,  these  derived  data  do  not  provide  additional  Information  about  the  system  dynamics; 
however,  they  help  to  keep  the  speed  response  of  the  model  within  a  realistic  range  and  to  prevent  long¬ 
term  speed  drifts.  This  characteristic  is  Important  since  the  identification  procedures  requires  the 
Integration  of  highly  unstable  (hover)  system  differential  equations  for  a  time  duration  of  about  90  sec. 

b.  Identification  of  the  Lateral/Directional  Motion 

Preliminary  time-domain  Identification  analyses  showed  that  the  longitudinal  and  lateral/directional 
motions  of  the  XV-15  are  practically  decoupled.  The  main  emphasis  was  placed  on  the  identification  of  a 
linear  lateral/directional  model.  This  model  is  represented  by  (linear)  differential  equations  for  the 
lateral  force,  rolling  moment,  and  yawing  moments.  The  general  3  DOF  model  is: 

x  -  Ax  +  Bu  +  bx 

(9) 

y  =  Cx  +  Du  +  by 


x  =  (v,  p,  r,  ♦) 
y|  ■  (a¥.  v,  p,  r,  ♦) 

uT  -  (s'.  sr) 

The  unknown  coefficients  in  the  state  matrices  (A  and  C)  and  the  control  matrices  (B  and  D)  are  the 
desired  stability  and  control  derivatives.  The  bias  vectors  bx  and  by  are  estimated  constants  repre¬ 
senting  drifts  and  zero  shifts. 

In  each  flight  test,  a  controller  (either  6a  or  Ar)  was  used  to  excite  the  aircraft  modes.  To 
obtain  sufficient  information  about  both  roll  and  yaw  motion  for  the  identification,  data  obtained  from  an 
aileron  and  a  rudder  sweep  were  combined.  This  "multiple  run  evaluation"  yields  one  common  model  for  both 
runs  (except  for  the  bias  terms,  which  must  be  estimated  for  each  individual  run).  This  approach  has  been 
used  successfully  In  previous  helicopter  identification  studies  (Ref.  5). 

Three  main  characteristics  of  the  XV-15  lateral/directional  dynamics  became  obvious  from  the  Initial 
identification  analysis: 


The  yaw  motion  which  is  due  to  rudder  inputs  is  virtually  decoupled  and  the  significant  parameters, 
yaw  damping  and  the  control  derivative,  can  easily  be  extracted  from  the  rudder-sweep  data.  Yaw- 
model  and  aircraft  time  histories  are  in  good  agreement. 

There  is  some  coupling  from  the  aileron  inputs  to  the  yaw  motion.  Therefore,  the  control -coupl ing 

derivative  N6  was  Included  for  identification. 

a 

For  the  aileron-sweep  data,  it  was  not  possible  to  obtain  a  satisfactory  curve  fit  for  the  total  run 
duration.  The  major  difficulty  is  the  identification  of  the  roll-moment  equation  and,  consequently, 
the  fit  of  the  roll -rate  response. 

The  third  characteristic  caused  some  severe  Identification  problems  and  was  investigated  in  more  detail. 

One  approach  to  this  problem  was  to  use  shorter  time  intervals  of  the  aileron  sweep  (only  the  low-  or 
mid-  or  high-frequency  part).  With  this  approach,  the  responses  of  the  identified  models  fit  the  measured 
roll  rates  almost  perfectly.  However,  there  were  major  differences  in  the  estimated  parameters  from  the 
original  Identification  based  on  the  total  run  duration.  Tests  wfth  different  a  priori  values  to  start 
the  ML  technique  were  also  made  to  ensure  that  the  ML  criterion  did  not  lead  to  local  minima  (a  conrion 
identification  problem).  Results  from  these  calculations  clearly  showed  that  the  data  contain  strong 
nonlinearities  which  cannot  be  described  by  a  linearized  model.  One  logical  next  step  is  the  extension  of 
the  model  to  include  the  appropriate  nonlinearities;  this  extension  will  be  addressed  later.  Another 
possibility  is  to  stay  with  a  linear  model,  accept  Its  deficiencies,  and  define  its  range  of  validity  and 
applicability.  This  approach  is  discussed  first. 

Lateral/directional  model  Identification  was  conducted  separately  on  the  three  available  aileron- 
sweep  runs,  each  in  combination  with  a  rudder-sweep  run.  When  the  total  run  duration  was  used,  all  three 
sweep  results  showed  the  same  tendency: 

The  model  response  matched  the  low-frequency  part  of  the  data  fairly  well. 

The  model  response  underestimated  the  flight  data  as  the  input  frequency  increased,  with  up  to  a  50% 
error  in  roll  rate  at  high-input  frequencies. 

These  results  make  sense  in  light  of  the  ML  identification  criterion: 

N  T  . 

I  -  E  ((2(1,)  -  y(t, ))  •  R'1  •  (z (t , )  -  y (t , ) ) I  ♦  m  ■  In  |R|  (10) 

where 

N  =  number  of  data  points 
z  =  measurement  vector 
y  =  model  response  vector 
R  =  measurement  noise  covariance  matrix 

The  optimum  is  reached  when  the  differences  between  the  amplitudes  of  the  measured  and  calculated  time 
histories  are  minimized.  From  fig.  4b,  it  is  seen  that  about  70%  of  the  total  run  duration  of  the 
aileron-sweep  is  low-frequency  data.  Consequently,  the  identification  method  emphasizes  primarily  the 
longer-duration,  low-frequency  part  of  the  data,  and  sacrifices  the  accuracy  of  the  shorter,  high- 
frequency  part.  For  many  applications,  the  Initial  and  short  term  (higher-frequency)  response  of  a  system 
is  of  more  Interest  than  the  long  term  (lower-frequency)  behavior.  Therefore,  It  was  desirable  to  improve 
the  identification  result  .'or  the  higher-frequency  range,  allowing  larger  errors  for  the  low  frequen¬ 
cies.  Methods  to  meet  this  objective  are: 

i.  Conduct  f requercy-sweeps  with  mor»  emphasis  on  the  high-f requenry  content. 

ii.  Apply  alternate  control  :nputs  (e.g.,  mu>t isteps)  which  excite  mostly  the  mid-  and  high- 
frequency  dynamics. 

iii.  use  only  the  higher -f requency  sweep  data  for  the  identification. 

'hese  approaches  were  either  not  possib  e  (new  'light  testing  reouired  for  options  1  and  ii)  or  they  were 
felt  to  be  a  poor  compromise  (Hi). 

Another  solution  is  to  ‘n'rease  tre  influence  of  tne  amplitude  errors  for  a  selected  part  of  tne 
data.  When  * requency - sweep  inpjts  are  jseu,  this  ran  be  lone  by  the  “multiple  segment  evaluation":  a 
part  of  the  data  (».g.,  hlgh-f -equency  range)  Is  treated  as  a  separate  test.  It  Is  combined  several  times 
with  the  original  test  data  so  that,  in  principal,  the  weight  ng  of  the  rhosen  data  points  is  arbitrarily 
increased.  This  approach  worxed  satisfactorily,  bet  it  y'e’ded  an  increased  number  of  unknown  biases, 
needed  more  data  handling  and,  in  particular,  required  more  computing  time.  But  pursu’ng  tr  y  basic  idea, 
the  identification  program  was  modified  to  allow  different  weighting  of  selected  time  periods  within  one 
run.  This  approach  turned  cut  to  b<  very  efficient  as  it  does  not  require  estimating  any  additional 
parameters,  or  computing  capacity. 
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The  data  weighting  technique  was  app’ied  tor  the  identification  of  the  aileron  f requency - sweeps . 
Increased  weights  were  used  for  the  roll-rate  fit  errors  which  occurred  in  the  nigner-frequency  part  of 
the  data.  Good  agreement  of  the  measured  and  calculated  data  was  thus  obtained  for  the  mid-  and  high- 
frequency  inputs,  whereas  there  were  larger  discrepancies  for  the  'ow- frequency  inputs.  Also,  there  was 
good  consistency  of  results  for  the  three  repeat  runs.  The  results  also  confirmed  that  't  's  advantageous 
to  keep  the  low-frequency  data  in  the  evaluation  as  they  provide  tne  necessary  speed -oer 1  vat i ve  ’nforma 
1 1  on .  The  mean  values  and  standard  deviations  for  the  derivatives  obtained  from  the  identification  of  the 
aileron  sweeps  are  summarised  in  Table  ?.  As  time-domain  techniques  tend  to  be  sens’t've  to  phase  sn'ts, 
a  time  lag  for  the  control  Input  was  estimated  as  a  multiple  of  tne  samp  i  ng  •■He.  in  state  space  format, 
the  final  time-domain  Identification  (mean  value)  mode’  for  tne  latera  i /qi  reel  'ona 1  mot’or  •%■_ 


units:  v  :  m/sec 

p.r  :  rad/sec 
*  :  rad 


V‘r  :  de9 

*Time  delay  in  control  input  is  *  0.03?C  sec. 

Figure  8  gives  the  time-history  comparison  for  one  of  tne  sweeps  w'f  tie  state  spa  r  v 
this  final  time-doma'n  model  correlates  well  at  med'um  and  h-gh  .-f  some  t  s  -e;a' 

frequency.  Overall,  however,  the  agreement  is  quite  sat ’ sf ac t or , . 


- flight  data 

- T IMF  DOMAIN  IDENTIFICATION 

10  r - 


TIME.  *er 

Fig.  8.  Time-domain  identification  of  1 atenal /d irect 1 onal  model. 
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The  present  XV- 15  study  Is  the  first  experience  with  explicit  data  weighting.  This  technique  cer¬ 
tainly  requires  further  development  and,  as  with  all  such  weighting  methods,  should  be  used  very  care¬ 
fully.  In  this  regard,  the  comparison  with  the  frequency-domain  results  Is  very  helpful  In  evaluating  the 
confidence  and  the  range  of  validity  of  the  results. 

c.  Identification  of  the  Longitudinal  Motion 

For  the  Identification  of  the  longitudinal  model,  an  elevator  sweep  test  was  combined  with  a  power 
sweep  test.  The  identification  results  for  the  longitudinal  dynamics  were  analogous  to  the  preceding 
lateral/directional  results: 

The  heave  equation  is  practically  decoupled  and  can  easily  be  identified  from  the  power  sweep  tests 
to  obtain  the  vertical  damping  and  the  control  derivatives. 

The  elevator  sweep  showed  the  same  tendency  as  the  aileron  sweeps:  It  was  not  possible  to  determine 
a  model  that  Is  equally  good  for  the  low-  and  high-frequency  range.  Again,  the  main  problem  occurred 
in  the  moment  equation  so  that  the  discrepancies  were  seen  In  the  pitch-rate  comparison. 

Only  one  of  the  three  available  flight  tests  could  be  evaluated.  When  the  other  two  tests  were  used,  the 
Identification  results  diverged  and  became  unusable.  This  Is  In  agreement  with  the  frequency-domain  anal¬ 
ysis  which  also  indicated  some  problems  with  these  runs.  For  the  one  remaining  elevator  sweep  run,  the 
first  30  sec  of  data  had  to  be  removed  In  order  to  reach  convergence  in  the  estimation.  Again,  the  data- 
weighting  technique  was  successfully  used  to  obtain  a  satisfactory  fit  for  the  higher  frequency  part  of 
the  data.  The  longitudinal  model  parameters  are  given  in  Table  2.  Unfortunately,  the  Identified  longitu¬ 
dinal  model  is  based  on  only  a  rather  limited  amount  of  data.  Therefore,  except  for  the  heave  equation, 
this  model  cannot  be  expected  to  have  the  same  level  of  reliability  as  the  lateral/directional  model. 
However,  the  good  comparison  with  the  frequency-domain  results  as  discussed  In  the  next  sections  show  that 
the  time-domain  model  accurately  represents  the  XV-15  longitudinal  dynamics. 


TABLE  2  Time-domain  Identification  Results 


Derivative  Mean  Variance  Standard  deviation  Standard  Deviation 

(%  of  mean  value) 


a)  Lateral/Directional  Parameters  (3  runs) 
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b)  Longitudinal  Parameters  (1  run) 
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4.  COMPARISON  Of  IDENTIFICATION  RESULTS 


This  section  compares  the  frequency-  ana  time-domain  identification  results.  This  comparison  is  done 
In  frequency-  and  time-domain  formats,  since  Doth  are  Important  for  ensuring  model  fidelity.  In  the 
frequency-domain  format,  transfer-function  parameters  from  the  frequency -domain  identification  are  com¬ 
pared  with  those  obtained  from  the  state-space  formulation.  Also,  frequency  responses  from  the  two  models 
are  compared  with  the  flight-data  frequency  response.  Since  the  freouency -domain  format  Is  the  "natural 
environment"  for  frequency-domai n  identification,  the  models  obtained  from  this  approach  generally  fit  the 
flight -data  responses  better  than  those  obtained  from  time-domain  identification.  Comparison  in  the  time- 
domain  format  is  achieved  by  driving  the  frequency-domain  models  with  the  frequency -sweep  input  histories. 
The  resulting  responses  are  compared  with  the  responses  of  the  vehicle  and  the  time-domain  identification 
fits.  Since  this  is  the  "natural  environment"  for  time-domain  ident if icat ion,  models  identified  using 
this  approach  generally  match  the  flight  data  better  here.  (The  detailed  discussion  of  the  results  for 
the  roll-axis  is  continued,  and  the  results  for  the  remaining  axes  are  again  summarized.) 

A.  Comparison  in  the  Frequency-Domain  Format 

Transfer  functions  are  obtained  from  the  time-domain  identification  results  of  Eqn.  11  (and  Table  2) 
by  Cramer's  rule,  and  are  tabulated  for  comparison  with  the  frequency-domain  results  in  Table  1. 

(1)  Lateral /directional  Models 

The  results  of  Table  1  indicate  that  the  lateral/directional  modes  (denominator  factors  of  the  trans¬ 
fer  functions)  are  nearly  identical  for  both  techniques,  except  for  the  difference  in  the  unstable  damping 
ratio  (cr).  The  high-frequency  gain  and  time  delay  of  the  three  transfer-functions  compare  very  well, 
while  there  are  some  differences  in  the  low-frequency  numerator  factors. 

The  relative  significance  of  the  differences  in  the  transfer-function  parameters  can  be  more  clearly 
seen  from  the  frequency-response  comparison  of  the  models  with  the  flight  data.  The  roll  responses  of  the 
identified  models  and  the  aircraft  are  shown  in  Fig.  9.  At  frequencies  of  u  >  1  rad/sec,  both  models 
correspond  almost  exactly  with  the  flight  data.  Also,  both  models  correctly  predict  a  low-frequency 
instability  at  w  *  0.5  rad/sec,  with  a  falling  magnitude  response  for  lower  frequencies. 

A  closer  examination  of  the  magnitude  and  phase  comparisons  shows  that  the  frequency-domain  identifi¬ 
cation  result  matches  the  magnitude-response  curve  better  than  the  time-domain  identification  result  in 
the  vicinity  of  the  dominant  mode  (u  *  0.5  rad/sec).  However,  the  time-domain  identification  result 
matches  the  phase-response  curve  better  in  this  frequency-range.  This  difference  is  due  entirely  to  the 
inherent  weighting  of  the  two  methods.  In  the  frequency-domain  identification  method,  the  relative 
weighting  between  magnitude  and  phase  is  arbitrary,  but  the  standard  choice  (I  DB  error:  7°  error)  has 
produced  satisfactory  results  In  a  number  of  identification  studies  conducted  by  one  of  the  authors 
(Refs.  1-4),  In  time-domain  Identification,  the  performance  index  is  much  more  sensitive  to  phase  errors, 
which  generate  a  large  area  between  the  model  and  flight -data  responses.  Thus,  the  phase  response  is  more 
closely  matched.  Also,  time-domain  identification  results  can  be  highly  sensitive  to  the  identified  value 
of  time-delay,  which  must  be  an  Integral  multiple  of  the  sample  rate. 

The  fact  that  magnitude  and  phase  curves  cannot  be  matched  simultaneously  in  either  frequency  or 

time-domain  identification  methods  further  Indicates  the  existence  of  important  nonlinearities  in  the  low- 

frequency  roll  oscillation  modes.  Therefore,  linear  models  (from  either  method)  are  a  compromise  and 

cannot  fully  characterize  the  nonlinear  behavior  of  the  vehicle.  Both  methods  capture  the  important 

vehicle  response  characteristics  and  are  generally  in  good  agreement  with  each  other.  Similar  agreement 

is  also  exhibited  in  the  yaw  responses  to  rudder  (r/i  )  and  aileron  (r/6  1  inputs. 
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Fig.  9.  Comparison  of  roll-response  models  (p/-ss)  in  the  frequency-domain  format,  a)  Magnitude; 
b)  phase. 


(2)  Longitudinal  Models 

The  comparison  of  transfer-function  models  for  the  longitudinal  degrees-of-freedom  is  very  similar  to 
the  preceding  results  for  the  lateral/directional  degrees-of -freedom.  The  dominant  modes  of  motion  for 
the  two  methods  are  very  close,  except  for  the  difference  in  the  unstable  damping  ratio,  tp  (again  roughly 
a  factor  of  2).  The  high-frequency  gain  and  time  delays  of  the  two  transfer  functions  are  also  nearly 
identical,  with  some  differences  in  the  low-frequency  numerator  parameters.  As  before,  the  frequency- 
response  match  between  the  two  models  and  the  flight  data  Is  nearly  identical  for  frequencies  greater  than 
l  rad/sec.  Also,  in  the  frequency  range  near  the  dominant  mode  (u  =  0.5  rad/sec),  the  frequency-domain 
model  fits  the  magnitude  response  better,  while  the  time-domain  model  fits  the  phase  response  better. 

Once  again,  nonl Inearities  and  differences  in  Inherent  weighting  of  the  methods  is  the  cause  of  this 
discrepancy.  In  general,  however,  the  agreement  between  the  models  and  the  flight  data  is  quite 
satisfactory. 

B.  Comparison  in  the  Time-Domain  Format 

Transfer  functions  obtained  from  the  frequency-domain  identification  were  converted  into  a  canonical 
state-space  representation  to  generate  time  histories  for  the  comparison  with  the  flight-test  data.  A 
bias  term  was  estimated  for  each  equation  (using  a  least-squares  procedure)  to  compensate  for  zero  shifts 
and  drifts.  Figure  10  shows  an  aileron-sweep  time  history  compared  with  the  frequency-domain  identified 
model,  and  the  time-domain  identified  model.  For  the  high-frequency  Inputs,  both  models  yield  virtually 
the  same  result  and  agree  with  the  flight-test  data.  In  the  lower  frequency  range,  some  differences 
between  the  two  models  and  differences  with  the  flight-test  data  can  be  seen.  Generally,  the  agreement 
with  the  flight  data  is  quite  satisfactory,  so  It  can  be  stated  that  both  Identified  models  represent  the 
dynamics  of  the  aircraft  fairly  well.  The  discrepancies  between  the  two  model  responses,  however.  Indi¬ 
cate  that  no  unique  model  can  be  identified;  the  slightly  different  results  reflect  the  specific  Identifi¬ 
cation  criterion  of  each  method.  This  confirms  the  preceding  conclusions  from  the  comparison  In  the 
frequency-domain  furmat. 


5.  TIME-00MAIN  VERIFICATION  USING  STEP-RESPONSE  DATA 

A  good  way  to  judge  the  utility  of  the  identification  results  Is  to  compare  the  prediction  of  the 
identified  models  with  the  vehicle  response  for  inputs  other  than  those  which  were  used  in  the  Identifica¬ 
tion  procedure.  Here,  step  Inputs  are  used  since  these  are  quite  different  from  the  frequency-sweep  forms 
which  are  used  in  Identification.  (These  step  inputs  tended  to  be  very  rounded  In  nature,  so  low-pass, 
preconditioning  to  remove  hlgh-f requency  elements  of  the  input  signal  Is  not  necessary  as  was  done  In 
Ref.  A.) 
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Fig.  10.  Comparison  of  ml ' -response  models  in  the  time-domain  fonnat. 


A.  Lateral/DIrectlona)  Motion 


The  Identification  of  the  yaw  motion  did  not  cause  any  difficulties.  Therefore,  a  satisfactory  pre¬ 
diction  capability  can  be  expected.  Figure  11  compares  the  yaw-rate  model  responses  with  the  measured 
data  for  a  pedal-step  input.  Good  agreement  is  apparent  for  both  models.  The  small  discrepancies  are 
caused  by  an  inaccurate  calibration  factor  between  pedal  and  rudder  (surface)  deflection,  and  some  mid¬ 
frequency  mismatch  of  the  first-order  yaw  model  (Ref.  9). 

Identification  problems  associated  with  the  aileron-sweep  evaluation  have  been  discussed  in  detail. 
The  verification  using  step  input  data  offers  a  good  possibility  to  check  the  validity  of  the  linear 
models.  Lateral  stick  step-inputs  were  flown  with  the  roll  SCAS-off  and  the  yaw  SCAS-on;  so  the  measured 
yaw-rate  response  shown  in  Fig.  12  results  from  the  pilot’s  lateral  stick  and  pedal  inputs,  and  the  yaw- 
SCAS  activity.  The  comparison  of  the  roll  rate  (p)  and  roll  angle  (*)  response  proves  that  both  models 
are  able  to  predict  accurately  the  aircraft  motion.  This  agreement  is  also  true  for  the  yaw  rate  compari¬ 
son  which  Indicates  that  the  coupling  derivative  (Nsa)  was  correctly  Identified.  Minor  differences 
between  the  two  model  responses  probably  result  as  before  from  the  different  weighting  methods. 

B.  Longitudinal  Motion 

The  heave  response  is  practically  decoupled  and  gave  no  problems  in  either  identification  method; 
good  verification  results  are  expected.  Figure  13  shows  that  the  power  step  responses  agree  with  the 
measured  (quite  noisy)  vertical  acceleration  data.  The  responses  are  shown  separately  for  the  two  models 
since  they  are  practically  Identical  and  cannot  be  distinguished  when  shown  within  the  same  plot. 
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Fig.  11,  Comparison  of  yaw-rate  response 
prediction  for  step  rudder  input  (yaw  SCAS-off). 
a)  Pedal  input  +  100%  spe<j  *  +44  deg  sr;  b)  yaw- 
rate  response. 


Fig.  12.  Comparison  of  lateral/directional 
response  prediction  for  step  aileron  input  (roll 
SCAS-off,  yaw  SCAS-on).  a)  Aileron  input;  b)  roll 
rate;  c)  roll  angle;  d)  yaw  rate. 
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Fig.  13.  Comparison  of  vertical-acceleration  response  prediction  for  step  power  lever  input,  a)  Power 
lever  input;  b)  frequency-domain  model;  c)  time-domain  model. 


Non linearities  and  a  limited  pitch-response  data  base  led  to  problems  in  identifying  a  longitudinal 
model,  as  has  been  discussed.  Therefore,  the  verification  tests  using  elevator  steps  are  particularly 
helpful  in  checking  the  reliability  of  the  two  models.  Figures  14  and  15  compare  the  pitch-model 
responses  for  two  different  flight  conditions:  pitch  SCAS-off  and  pitch  SCAS-on.  In  both  cases,  the 
identified  models  yield  a  good  prediction  of  the  aircraft  response.  Again,  the  minor  differences  between 
the  responses  are  due  to  the  inherent  weighting  of  each  method. 

The  SCAS-off  data  fit  of  Fig.  14  is  of  special  interest.  Since  the  models  were  extracted  from 
SCAS-on  flight-test  data,  some  output/input  correlation  cannot  be  avoided  and  may  lead  to  significant 
identification  errors;  in  the  worst  case,  the  inverse  feedback  transfer  function  rather  than  the  open-loop 
aircraft  response  would  be  identified  (Ref.  9).  However,  the  good  agreement  between  the  model  time  his¬ 
tories  and  the  SCAS-off  flight  data  in  Fig.  14  clearly  demonstrates  that  the  open- loop  dynamics  of  the 
aircraft  were  determined.  The  overall  excellent  correlation  of  the  models  and  step-response  data  adds 
confidence  to  the  accuracy  of  identified  derivatives  and  transfer  functions,  and  the  estimation 
techniques. 


6.  NONLINEAR  MOOEl  IDENTIFICATION 

The  preceding  identification  results  from  both  time-  and  frequency-domain  techniques  have  demon¬ 
strated  that  linear  model  Identification  yielded  a  compromise  between  low-  and  high-frequency  data  fits, 
or  between  magnitude-  and  phase-response  fits.  They  suggest  the  existence  of  significant  nonllnearltles, 
particularly  in  the  roll  and  pitch  axes.  Therefore,  the  nonlinear  maximum  likelihood  time-domain  method 
was  utilized  to  identify  an  extended  model  and  to  investigate  the  importance  of  various  parametric  terms. 

Relatively  large  amplitude  aircraft  responses  during  the  low-frequency  Inputs  (see,  for  example, 

Fig.  8)  are  a  common  characteristic  of  all  of  the  frequency- sweep  flight-test  data.  Deviations  from  the 
steady-state  trim  are  in  the  range  of: 

9-14  m/sec  In  lateral  speed  (aileron-sweep) 

8-11  m/sec  in  longitudinal  speed  (elevator- sweep) 

25-37  deg  in  roll  angle  (aileron-sweep) 

17  deg  in  pitch  angle  (elevator-sweep) 
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Fig.  14.  Comparison  of  pitch-response  prediction 
for  step  elevator  Input  (pitch  SCAS-off). 
a)  Elevator  input;  b)  pitch  rate;  c)  pitch  angle. 
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Fig.  15.  Comparison  of  pitch-response  prediction 
for  step  elevator  Input  (pitch  SCAS-on). 
a)  Elevator  Input;  b)  pitch  rate;  c)  pitch  angle. 


These  amplitudes  are  certainly  large  enough  to  violate  the  small-perturbation  assumptions  for  linear 
models;  further,  the  dynamic  characteristics  of  hovering  rotorcraft  are  especially  sensitive  to  transla¬ 
tional  speed  changes.  Therefore,  the  linear  model  was  first  extended  by  adding  nonlinear  speed  deriva¬ 
tives  (l(v**2)  and  L(v**3))  to  the  roll  moment  equation:  the  curve-fits  Improved  significantly.  In  par¬ 
ticular  owing  to  L(v**2).  Further  attempts  to  reduce  the  remaining  discrepancies  were  made  by  Including 
additional  nonlinear  terms.  Their  significance  was  checked  with  time-history  comparisons  and  the  evalua¬ 
tion  of  the  parameter  covariance  matrix,  which  Indicates  the  reliability  of  the  identified  parameter  and 
the  correlation  with  other  parameters.  As  a  preliminary  result,  a  model  was  Identified  that  Includes  the 
above  mentioned  speed  derivatives  and,  In  addition,  L(a#**2)  and  L(6a*u). 

Figure  16  shows  that  the  nonlinear  model  fits  the  measured  data  almost  perfectly.  The  results  pre¬ 
sented  In  Fig.  16  are  preliminary  and  are  Intended  to  Illustrate  the  possible  role  of  nonl inearltles  in 
the  dynamics.  It  Is  Important  to  note  that  the  model  was  Identified  without  the  use  of  any  explicit  data 
weighting.  This  suggests  that  the  additional  weighting  (e.g.,  high  frequency  versus  low  frequency)  Is  not 
required  when  an  appropriate  model  formulation  Is  applied.  However,  the  evaluation  again  revealed  a  well 
known  identification  problem:  It  Is  always  possible  to  Improve  the  time-history  curve  fit  by  arbitrarily 
adding  model  parameters.  But,  a  useful  model  requires  the  estimated  derivatives  to  have  physical  signifi¬ 
cance.  It  Is  the  responsibility  of  the  analyst  to  define  and  select  meaningful  additional  terms.  For  the 
side-by-side  rotor  configuration  of  the  XV-I5,  the  speed-related  derivatives  (l(v**2)  and  L(v»*3))  are 
physically  Justified.  Similarly,  the  control  effectiveness  ma^  be  In  fact  nonlinear  and  dependent  on 
forward  speed;  but  these  effects  should  be  further  Investigated. 


7.  ASSESSMENT  OF  IDENTIFICATION  METHOOS 

This  cooperative  study  has  provided  the  unique  opportunity  for  specialists  using  different  methods  to 
compare  and  coordinate  analyses  of  a  common  rotorcraft  data  base.  This  experience  has  been  Invaluable  for 
gaining  a  better  appreciation  for  the  advantages  and  limitations  of  both  techniques,  and  for  formulating 
Ideas  for  an  Integrated  approach  to  dynamics  Identification.  The  following  observations  are  based  on  the 
results  of  this  cooperative  effort. 
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Fig.  16.  Example  of  time-domain  Identification  of  a  nonlinear  lateral/dlrectlonal  model. 


The  principal  advantages  of  frequency-domain  Identification  are: 

a.  Frequency  responses  of  the  dominant  on-axis  Input/output  pairs  are  rapidly  generated  and  are  very 
useful  for  gaining  a  good  appreciation  for  the  Inherent  vehicle  dynamics.  The  fact  that  the  extracted 
frequency-responses  are  Independent  of  pre-assumed  models  is  Important  for  the  Initial  assessment  of 
natural  system  orde r,  dominant-mode  locations,  and  stability  characteristics.  As  a  result,  a  better 
choice  of  appropriate  model  structure  and  order  is  possible. 

b.  Parameters  associated  with  the  high-frequency  dynamic  behavior  can  be  determined  directly  from 

the  frequency-responses  without  any  a  priori  assumption  of  model  structure.  Specifically,  the  control 

derivatives  (e.g.,  L .M*  )  are  determined  from  the  high-frequency  gain  responses,  and  the  equivalent  time 

delays  (e.g.,  t  ,t  )  are  determined  from  the  high-frequency  phase  responses. 

•  e 

c.  Weighting  can  be  accomplished  explicitly.  Relative  weights  can  be  arbitrarily  assigned  to  the 
magnitude  and  phase  curves.  Model  fitting  can  also  be  arbitrarily  weighted  more  to  the  low-  or  high- 
frequency  range- -depending  on  the  Intended  use  of  the  model. 

d.  Accurate,  hlgh-resolutlon  frequency-response  Identification  is  given  the  main  emphasis  in  this 
method.  CHIRP  ?-transform  methods  are  especially  well  suited  for  Identifying  f requency-responses  from 
noisy  flight  data.  The  resulting  transfer-function  models  are  a  much  closer  representation  of  the 
frequency-response  characteristics  than  Is  possible  with  time-domain  Identification. 

The  principal  disadvantages  of  frequency-domain  Identification  are: 

a.  Current  techniques  are  not  well  suited  for  highly  coupled  multi- tnput/mul tl -output  (MIMO)  system 
identification,  although  two-lnput/slngle-output  Identification  has  been  successfully  attained  In  the 
present  study.  More  highly  automated  techniques  are  needed  to  make  the  frequency -domain  methodology  effi 
clent  when  the  required  number  of  Input/output  frequency-responses  Is  large.  Also,  methods  for  slmulta 
neous  fitting  of  many  coupled  responses  is  necessary  to  ensure  commonality  of  transfer-function  denomi¬ 
nator  parameters  for  MIMO  models. 


b.  Frequency -domain  identification  results  in  transfer-function  models.  Individual  stability  deriv¬ 
atives  are  not  readily  extracted  unless  the  assumed  models  are  of  very  low  order. 

c.  Spectral  analysis  assumes  input -to-output  linearity.  For  nonlinear  systems,  the  transfer  func¬ 
tions  are  linearised  describing  functions;  identification  of  pure  nonlinear  parameters  is  not  possible. 

The  principal  advantages  of  time-domain  identification  are: 

a.  The  method  is  naturally  suited  to  multi-input/multi-output  identification  since  the  model  can  be 
of  arbitrary  order  and  structure.  This  method  is  especially  well  suited  for  Identifying  highly  coupled 
systems. 


b.  Stability  and  control  derivatives  are  identified  explicitly,  and  the  method  leads  to  the  identi¬ 
fication  of  a  complete  state-space  model. 

c.  Considerable  effort  is  invested  in  achieving  the  highest  quality  of  time-domain  data.  Data  con¬ 
sistency,  drop-out  tests  and  signal-reconstruction  methods  are  an  Integral  part  of  the  time-domain  identi¬ 
fication  procedure.  The  least-squares  fitting  in  the  time-domain  with  high-quality  time-history  data 
results  in  a  much  better  time-domain  fit  of  the  frequency- sweep  responses. 

d.  Extended  maximum- 1  ike  1 i hood  techniques  can  be  used  to  identify  parametric  nonlinearities  which 
are  especially  Important  in  the  low-frequency  dynamics  of  hovering  rotorcraft. 

The  principal  disadvantages  of  time-domain  Identification  are: 

a.  The  results  are  dependent  on  presumed  model  structure  and  order,  when  a  new  configuration  is 
being  identified,  a  priori  knowledge  of  model  structure  may  not  be  available,  and  considerable  variation 
in  the  parameters  can  occur  when  the  model  structure  is  altered. 

b.  Explicit  frequency-domain  weighting  is  not  possible.  Specifically,  the  time-domain  method  inher¬ 
ently  weights  phase  errors  more  heavily  than  magnitude  errors.  This  characteristic  makes  the  extracted 
state-space  model  very  sensitive  to  pure  time  delays  and  unmodeled  high-frequency  dynamics.  Also,  the 
method  inherently  weights  low-frequency  dynamics  much  greater  than  high-frequency  dynamics;  this  weighting 
can  be  adjusted  when  the  input  signal  has  monotonical ly  increasing  frequency  content  as  in  the  frequency- 
sweep. 


c.  Confidence  in  the  individual  state-space  model  parameters  may  be  very  low  since  the  identified 
st»*e-space  model  can  contain  a  high  degree  of  Internal  cancellation  in  the  overall  Input-to-output 
response. 


8.  A  PROPOSAL  FOR  A  COORDINATED  FREQUENCY -D0HA1N/T1ME -DOMAIN  IDENTIFICATION  METHOD 

The  preceding  assessment  of  the  advantages  and  limitations  of  each  identification  method  suggests  the 
following  coordinated  frequency-domain/tlme-domain  Identification  approach: 

Step  1.  Use  time-domain  signal  conditioning  methods  to  clean  up  the  flight  data  for  drop-outs,  wild 
points,  and  consistency.  For  example,  rate  gyros  can  be  Integrated  and  compared  with  attitude  gyros. 

Step  2.  Identify  the  dominant  Input/output  on-axis  frequency-response  characteristics  using  only  the 
Dest  runs,  as  determined  from  coherence  analyses  of  the  individual  frequency-sweeps.  Identify  the  effec¬ 
tive  time  delay  and  high-frequency  control  sensitivity  directly  from  the  frequency-response  plots. 

Step  3.  Formulate  low-order  system  models  from  inspection  of  the  identified  frequency-response  plots 
and  theoretical  analyses.  Determine  the  on-axis  transfer-function  parameters. 

Step  A.  Formulate  a  state-space  model  which  has  a  structure  and  order  consistent  with  the  transfer- 
function  model  formulation.  Time-domain  identification  should  be  completed  with  the  equivalent  time  delay 
fixed  at  the  value  identified  in  Step  3.  Weighting  should  be  applied  to  the  time-history  data  to  ensure 
that  the  control  sensitivity  derivatives  are  maintained  at  the  value  identified  in  Step  3.  (Alterna¬ 
tively,  the  control  derivatives  can  be  fixed.) 

Step  5.  Compare  the  extracted  on-axis  transfer-functions,  frequency-responses  and  time  histories 
from  the  time-domain  and  frequency -domain  results.  If  substantial  errors  exist,  go  back  to  Step  2; 
reevaluate  the  quality  of  the  spectral -responses,  time  responses,  and  the  order  and  structure  of  the 
selected  models.  If  the  models  are  found  to  be  the  best  which  can  be  achieved  under  the  assumption  of 
linearity,  pursue  nonlinear  maximum- 1  ike  1 i hood  methods  to  identify  the  dominant  parametric  nonlinearities. 

Step  6.  Verify  the  extracted  models  using  time-history  data  from  inputs  not  used  in  the  identifica¬ 
tion  procedure.  If  significant  errors  between  the  predicted  and  actual  response  characteristics  exist, 
^evaluate  the  significance  of  observed  discrepancies  in  frequency  and  time-domain  identification  fits. 

If  necessary,  go  to  Step  3  and  increase  the  order  of  the  models;  but,  check  for  the  possibility  of  model 
over  parameterization  by  trying  a  few  different  verification  inputs. 
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9.  CONCLUSIONS 

This  joint  effort  has  provided  the  unique  opportunity  for  specialists  In  different  techniques  to 
compare  their  approaches  using  a  common  flight  test  data  base.  On  the  basis  of  this  comparison,  it  has 
been  shown  that  the  frequency  and  time-domain  methods  each  have  Important  advantages  and  Inherent  limita¬ 
tions.  Future  identification  efforts  must  be  based  on  a  comprehensive,  coordinated  approach  which  uses 
both  frequency  and  time-domain  methods. 
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